BNL-GP-1/93 
January, 1993 



Partition function for the eigenvalues of the 

Wilson line 



Andreas Gocksch 
& 

Robert D. Pisarski 

Department of Physics 
Brookhaven National Laboratory 
Upton, New York 11973 



Abstract 

In a gauge theory at nonzero temperature the eigenvalues of the Wil- 
son line form a set of gauge invariant observables. By constructing the 
corresponding partition function for the phases of these eigenvalues, we 
prove that the trivial vacuum, where the phases vanish, is a minimum of 
the free energy. 



In computing the properties of gauge theories at a nonzero temperature T, a stan- 
dard method of calculation is to use the imaginary time formalism, where the gauge 
fields are strictly periodic in a euclidean time r, with period (3 = 1/T JE]. This 
allows for a type of Aharonov-Bohm effect at nonzero temperature, characterized by 
a nonzero expectation value for the Wilson line that wraps around in the direction of 
the imaginary time, exp(ig J A Q dr). At one loop order [Q-Q] these inequivalent vacua 
are parametrized simply by a constant value for the gauge potential A ; calculation 
shows that the trivial vacuum, with Aq = 0, minimizes the free energy. At two loop 
order || it is necessary to chararcterize the vacua not just by constant Aq, but by 
the Wilson line itself after doing so, the trivial vacuum remains a minimum. 

As conjectured in ref. 0, in this paper we show that the trivial vacuum is always 
a minimum of the free energy. Our starting point is a gauge invariant quantity in 
euclidean spacetime, the effective potential for the phases of the eigenvalues of the 
Wilson line. This euclidean path integral is then transformed into a partition function 
||, where the constraint on the phases of the eigenvalues becomes an imaginary chem- 
ical potential for global color charge. Because this chemical potential is imaginary, 
it follows directly that the trival vacuum minimizes the free energy, up to standard 
degeneracies. 

Typically the partition function of a gauge theory is a sum not over all states, 
but only over those which are gauge invariant [|]||||. This is true, for example, in 
computing the free energy at nonzero temperature. In our case, while the sum over 
states is gauge invariant, the individual states which contribute are not. In this vein, 
we note that while partition functions with a chemical potential for global charge 
have been studied previously [|K|] , hitherto it was assumed that only gauge invariant 
states contribute. The extension of what is an allowed partition function is, we feel, 
the most striking feature of our results. 

For definiteness we consider an SU(N) gauge theory without matter fields, fol- 
lowing previous conventions [0. Under a local gauge transformation Q, the gauge 
potential transforms as A^ = W D M / \—ig) , with g the coupling constant, and 
Dfj, = <9 M — igA^ the covariant derivative. The Wilson line in the direction of eu- 
clidean time is 

L(x) = V exp(^ig A (x,t) oIt^J = U(x) A(x) U\x) ; (1) 
V denotes path ordering, and x is the coordinate for three spatial dimensions. The 
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matrix L(x) is unitary, and so we can write it as the unitary transformation of a 
diagonal matrix, A(x): Aij(x) = SijA^x), UW = 1. Under a gauge transformation 
the Wilson line transforms as 

L(x) -> Q(x, (3) V exp [ig A%(x, r) dr^j fi f (a;, 0) . (2) 

Without loss of generality we require that both the gauge fields and the gauge trans- 
formations are strictly periodic in the euclidean time. If Q(x,(3) = Q(x,0), then for 
the Wilson line a gauge transformation is just a similarity transformation; the matri- 
ces U(x) are gauge variant, U(x) — > Q(x,0)U(x), while the diagonal matrix A(x) is 
gauge invariant. 

Notice that what we call the Wilson line is a matrix in color space. This is is dis- 
tinct from what is usually termed Wilson (or Polyakov) loop at nonzero temperature, 
which is the trace of this matrix Jl| . Since tr L = tr A, the trace of the Wilson line is 
automatically gauge invariant. 

The phases of the eigenvalues are given by A(x) = e iX{x) ] the A are elements of the 
Cartan sub-algebra, which is the set of mutually commuting generators in the group. 
We introduce an effective potential for the phases as 

Z(~X) = I VA,(x,t) e- s s(\- f^X(x)) . (3) 

JA,,(x,l3)=+A„(xfi) \ J V ) 

VA^ is the functional measure for the gauge fields, including gauge fixing and ghosts, 
S is the gauge field action, and V is the volume of space. We have chosen to use a 
"constraint" effective potential |ll|; a more standard form, using an external source 



and then Legendre transformation, could also be used. From Z(X) the free energy, as 
a function of the A, is F(X) = —log(Z(X)). 

When the phases vanish, F(0) reduces to the usual free energy at a temperature 
T, in the absence of any background field. (This is seen most easily from ( jl^ ) below.) 
It is for this reason that we constrain the phases of the eigenvalues instead of the 
eigenvalues themselves. Since the A's are phases, the free energy is periodic: for 
Ajj = XiSij, F(X) is unchanged when any single element is shifted by a multiple of 2tt, 
F(A, + 2vr) = F(A l ). 

In a pure gauge theory there is a further degeneracy in F(X). The gauge poten- 
tials remain periodic under certain aperiodic gauge transformations, as long as the 
aperiodicity is a constant element in the center of the group For an SU(N) gauge 
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group the center is Z(N) , and the aperiodic gauge transformations satisfy Q(x,/3) = 
e £ Q(x, 0), where 0t = 2n£/N, and i is an integer = 0, 1, ...(N — 1). Under such gauge 
transformations every phase transforms by the same constant, Xi(x) — ► Aj(x) + 0t 
for all i = 1...N. This implies that the effective potential automatically posseses a 
Z(N) symmetry, F(X + 6t) = F(X) (here 8e denotes a diagonal matrix where each 
element = 9g; this notation is used later, and should be clear from the context). 

Before turning the effective potential into a partition function, it is worth review- 
ing the results of perturbative calculations To compute one expands about a 
constant, background field A = Aq + Aq U , where Aq = X/g. A convenient gauge is 
covariant background field gauge, where the gauge fixing term in the lagrangian is 



^covariant £ ( ("^A 1 ) ) ' 



D'rj = — ig[A c J:,}. At one loop order the free energy is independent of the gauge 
fixing parameter £, F oneloop (X) = 2tr log( — (D^) 2 ); this is the free energy for two 
spatially transverse gluons in the background field Aq. The minimum is at A = 0, 
F neioop{$) < F oneloop (X), plus Z(N) transforms thereof. 

At two loop order, as a function of Aq the free energy is gauge dependent. Belyaev 
0.0] noted that this happens because at one loop order fluctuations in the gauge field 
- in particular the static components of Aq U which do not commute with Aq — alter 
the relationship between A and Aq from that at tree level. This feeds back into the 
free energy at two loop order. As a function of the A, the free energy is independent 
of £, with a minimum for A = 0, F two i oop (0) < F twoloop (X). 

There are "unitary" gauges in which the relationship Aq = X/g is unchanged to 
any loop order. A reasonable guess is static gauge [|1^], imposing doA = 0. This is not 
adequate, however, for static gauge has a residual gauge freedom of performing static 
gauge transformations. Dividing up A = Aq+Aq U , some static gauge transformations 
rotate the background field Aq, which is inconsistent. A gauge which avoids this is 
static background gauge, obtained by adding the term 

Cstatic = ±trf(±D d >ir + Mr) ) (5) 

to the lagrangian. Letting the gauge parameter £' — >• fixes DqA^ = 0. This 
restricts all fluctuations in A^ to be both static and in directions which commute 
with Aq. The value of the other gauge fixing parameter £ is arbitrary; £ removes the 
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remaining gauge degeneracy under static gauge transformations which commute with 

Acl 
A Q . 

The calculation of the free energy to one loop order in static background gauge is 
illuminating. At leading order the ghost determinant depends on the background field 
as det(—(DQ 1 ) 2 ). In static background gauge the square root of the ghost determinant 
equals a product of Vandermonde determinants for the Wilson line, 

det(D$) = J] |V(A)| 2 , V(A) = II (e^-e^) . (6) 

The square of the Vandermonde determinant enters at each point in space, so in the 
action the logarithm of the product of Vandermonde determinants is proportional to 
5 3 (0)V, which vanishes in dimensional regularization. Even in other regularization 
schemes, however, the product of Vandermonde determinants from the static back- 
ground ghosts do not contribute to the free energy F(X). One factor cancels against 
the contribution of the spatially longitudinal mode of the gluon.[] The second fac- 
tor cancels against terms from the A components of the propagator, which give a 
finite but nonzero contribution in the limit £' — > 0. After these cancellations the two 
spatially transverse modes give the same result as found in covariant gauge. 

To turn the path integral of (|3|) into a partition function we follow Rossi and Testa 
||. First we transform to A = gauge. The gauge transformation which accom- 
plishes this transformation is the "partial" Wilson line, Q(x, r) = Vexp(ig J r A (x, r')dr'). 
Many factors can be ignored in going to Aq = gauge. The Fadeev-Popov deter- 
minant is det(d ), which is independent of temperature and the background field. 
Similarly, when < r < j3, all terms in the measure from the gauge transformation 
Q(x, r) can all dropped. The sole exception is the gauge transformation for the last 
slice of euclidean time, Q(x, f3). This must be retained because it affects the boundary 
conditions: the boundary conditions are aperiodic, Ai(x,/3) = +Af ^'(x, 0), while 

1 This cancellation is similar to that found by Weiss j2|, who imposed the condition that A$ is 
static and diagonal by fiat. There is no ghost per se, but the measure of the functional integral 
includes a single factor of the Vandermonde determinant squared at each point in space. He obtains 
the correct free energy at one loop order, where the product of Vandermonde determinants in 
the measure cancels against the contribution of the spatially longitudinal gluons. One can check, 
however, that his method does not give the correct form of Gauss' law in the Cartan sub-algebra. 
The problem is that the gauge condition 8qAq = can only be imposed after allowing arbitrary 
variations in A , not before. Imposing static background gauge as in the text, by letting £' — » 0, 
avoids this difficulty. 
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the Wilson line L(x) — > Q(x,/3). Relabeling Q(x,/3) as Q(x), in A = gauge the 
path integral is 

z(X) = /wk*) /^ MJfW) ^(».-) ^ >(i-f*m) • w 

As before we decompose = ?7(x)A(a;)C/^(a;), A(sc) = e tX< ^ x \ In terms of these 

variables the measure T>Vl = T>U T>\ |V(A)| 2 , where the measure naturally includes 
the square of the Vandermonde determinant Thus in Aq = gauge the constraint 
on the Wilson line becomes a constraint on the final gauge transformation, Q(x). 

An advantage of A = gauge is that the canonical structure is manifest 
Thus the functional integral in (^), written in imaginary time, is equal to the partition 
function, 



Z{\) = [ Vn(x) £ (A t (x)\ e-f™ \A?(x)) 5 ( A - / 

Ai(x) \ 



All fields Ai(x,t) are functions of a fixed time i; this time is arbitrary, and so the 
dependence upon it is suppressed. We assume canonical commutation relations be- 
tween Ai(x) and the electric field Ej(x), [Af(x), Ej(y)] = i 5 ab 5^ 5 3 (x — y). H is the 
total Hamiltonian, with the sum over all eigenstates \Ai(x)); \Af(x)) is the gauge 
transformation of a state under Q(x). 

To obtain a more useful form of the partition function we introduce the generators 
of gauge transformations. Given Q(x), let Q = e tuJ , and define the operator 

G(u) = -- I d i xtr((D, i pj(x))EAx)) . (9) 
g J 

The canonical commutation relations imply that these generators form a representa- 
tion of the Lie algebra, 

[GW.GN] = G{[lo u lo 2 )) . (10) 

By the operator identity Af(x) = e~ ^ Ai(x)e the gauge transformed state is 
given by 

\A?(x)) = e iG W\Mx)). (11) 

Since the operators G{u) generate gauge transformations, they commute with the 
Hamiltonian, [G(u),Tl] = 0. 
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We write the eigenvalues X(x) = A + X q (x), so that the constraint becomes 
/ d 3 xX q (x) = 0. This constraint implies that there is no constant mode in X q (x), 
and is automatically satisfied by requiring that X q (x) vanish at spatial infinity. In 
this way we trade the constraint for a boundary condition, 

Z(X) = I VU(x)VX q (x)\V(X+X q (x))\ 2 Yl ( A f( x )\ e^ n+iG ^ +iG ^ \Af{x)) , 

JX q (oo)=0 A . {x) 

(12) 

To obtain this we start with \A n (x)), and by using Q = UKU\ undo the last two 
gauge transformations. W is carried through to change (Ai(x)\ into (Af(x)\. The 
generators for the gauge transformations in the Cartan sub-algebra are left as is. The 
factor of G(X q (x)) generates a local, that for G(X) a global, gauge transformation in 
the Cartan sub-algebra. The generator of the global gauge transformation can be 
written as a chemical potential for global color charge, Q: 

G(X) = 2 tr(X Q) , Q = - ig jd 3 x[Mx),E t (x)]. (13) 

The appearance of the background field as an imaginary chemical potential is not 
surprising: this is obvious at tree level, using the canonical formalism to expand 
about a constant field Aq = X/g; see, e.g., ref.'s |§ and [j|]. What is not obvious, 



and which (|l^) and fll4]) demonstrate, is how to write the partition function Z(X) in 
a gauge invariant manner beyond tree level. 

It is worth checking that to leading order the partition function in ( |12"D gives the 
same free energy as found in the imaginary time formalism. The integrations over the 
gauge transformations U(x) and X q (x) can be ignored to this order. Thus the only 
change from calculating the free energy in zero field is the presence of an imaginary 
chemical potential for global color charge. Partition functions of this form have been 
studied from these calculations we see that the spatially transverse gluons give 
the same free energy as found to one loop order in the imaginary time formalism. The 
spatially longitudinal mode, with zero energy, cancels the product of Vandermonde 
determinants from the measure, l\ x \V(X + X q (x))\ 2 ~ H x |V(A)| 2 . 

Since X q vanishes at spatial infinity we can integrate G(X q ) by parts to obtain 
G(X q ) = 2/ d 3 x tr(X q DiEi)/ g . This presumes that the gauge field itself satisfies the 
appropriate boundary conditions: Ai(x) ~ 1/r and Ei(x) ~ 1/r 2 as x — > oo suffice. 
The only other place where X q enters is through the infinite product of Vandermonde 
determinants, Yi x + ^( x ))\ 2 - As remarked previously, however, this infinite prod- 
uct vanishes if we adopt dimensional regularization, which we do. Then integration 
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over X q is trivial, generating a projection operator, Pc = Yl x 5(tr(tcDiEi(x))) , where 
the are the N — 1 elements of the Cartan sub-algebra. This imposes Gauss' law in 
the Cartan sub-algebra at each point in space. Whence 

Z{\) = J VU(x) £ ( A i( x )\ e- pn + 2ttrCxQ) P c 1^0*0) • ( 14 ) 

Ai(x) 

This form of the partition function involves only states with positive norm, the 
\A^ x \x)}, and a positive measure of integration over U(x). The background field 
enters only through the global color charge, e 2 ** r ( A 2). While this factor is imaginary, 
because the sum always includes states with equal and opposite charge, the total is 
necessarily real, so e 2itr{XQ) ~ cos (2tr(A<2)). Since the cosine function is less than or 
equal to one, a nonzero field A decreases the partition function, Z(X) < Z(0), and so 
increases the free energy, 

F(0) < F(X) , (15) 

which concludes our proof. 

Our argument is inspired by that of Vafa and Witten [13] for the #-angle in QCD. 
Nonzero 9 appears in the euclidean path integral as a phase factor, e ieQtop , where 
Qtop is the topological charge. If the functional measure is positive definite, this 
implies that 6 = and tt minimize the corresponding free energy. In our case we 
need positivity for the sum over states in the partition function (which is why we 
take A = gauge) as the background field appears as a phase factor for global color 
charge. 

A delicate point is the use of dimensional regularization to eliminate the infinite 
product of Vandermonde determinants. Based upon calculations to one loop order, 
where the product of Vandermonde determinants cancels in -F(A), we suggest that 
this cancellation persists to arbitrary loop order. There is a general reason why it 
should. Consider Z(X) in (|l|), and then integrate over all A; the delta function drops 
out to give the free energy in zero field. (Integrating ( |T4"D with respect to A gives 
a constraint of zero color charge; since in infinite volume the minimum is always 
at A = 0, this constraint is inconsequential.) But the appropriate measure for the 
integration over A, V\ |V(A)| 2 , already includes a single factor of the square of the 
Vandermonde determinant. Thus there should not be any additional factors in Z(X), 
infinite or not. 

It is trivial to extend the above to include the effects of matter fields. Fermion 
fields, for instance, are added to the euclidean functional integral, (13) , through an- 
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tiperiodic boundary conditions. ip(x, j3) = —ip(x,0). The above derivations go 
through with minor modifications: the states include both gauge fields and fermions, 
\Ai(x),ip(x)), and the global color charge includes the contribution of the matter 
fields as well. But again we conclude that F(0) < F(X). 

The partition function also provides insight into the global Z(N) symmetry. If all 
matter fields lie in the adjoint representation, the free energy, and so the vacua, are 
Z(N) symmetric, F(0) = F(9 t ) for all £ = 0,1, ...(N - 1). This Z(N) symmetry is 
manifest from the form of the partition function in flU]). The factor of e 2ttr ^ deQ ^ = 
e iG{e e ) g enera t es a global Z(N) gauge transformation on the states by f2 = e ldl . As 
an element of the center of the group this transformation leaves gluon (and other 
adjoint) states unaltered, e 2itr ^\Af) = \Af). 

If there are matter fields in the fundamental representation, the Z(N) symmetry 
is lifted at one loop order: for £ ^ 0, F(0) < F(9i), as the trivial vacuum, A = 0, is 
the unique vacuum (modulo the usual factors of 2ti because A is a phase). In terms 
of the partition function, the Z(N) degeneracy is lifted because while states in the 
adjoint representation don't change under Z(N) global gauge transformations, states 
in the fundamental representation do: 

e 2 ' 1 * 2 '!^,/) = IAW> + \Af,^ u ). (16) 

In this way, we see that while the Wilson line is defined originally in euclidean space- 
time, its transformation under (aperiodic) global Z(N) gauge rotations directly re- 
flects the transformations of states, in the partition function, under global Z(N) gauge 
rotations. 

We conclude by discussing the nature of states which contribute to the partition 
function in (fH|). Since by definition Z(X) is gauge invariant, the sum in states in ( |H] ) 
must be as well. A (static) gauge transformation V(x) changes Ai — > A\\ both the 
bra and ket states change, (Af\ — > (Af u \ and \Af) — > |^4| /c/ ) . As the Haar measure is 
invariant under left multiplication, though, by U(x) — > V Jf (x)U(x) we can eliminate 
any such V(x). 

While the sum over states is gauge invariant, the states which contribute are not. 
This is very different from the computation of the free energy in zero field, F(0). 
Begin with (P, and just drop the constraint altogether. Then F(0) involves states of 
the form 

\A(x)} inv = Jvn(x) \A?(x)}. (17) 
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By the properties of the Haar measure, these states are manifestly gauge invariant 
IHH, \A\) inv = \Ai) inv . In nonzero field, however, the presence of the constraint on 
the eigenvalues implies that more than just gauge invariant states contribute to the 
sum. 

Alternately, consider how Gauss' law works. In zero field it is easy to start from 
(§) and show that the integral over Q imposes Gauss' law in all directions of the 
group This can be seen, albeit less directly, from (|T1|): for A = 0, take each 

state, say the bra (Af\, and undo the gauge transformation U(x) by writing it as 
{Ai\e~ tG ( u \ U = e m . Doing this for both bra and kets, the projector for the Cartan 
sub-algebra becomes e~ tG< ^~P ' c e tG( ^> '. This is the gauge transformation of Gauss' law 
in the Cartan sub-algebra, so integration over the U(x) imposes Gauss' law in all 
directions of the group algebra. 

These manipulations fail in nonzero field, A 7^ 0. The problem is that the generator 
for the local gauge transformation, G(u), need not commute with that for the global 
gauge transformation, G(X): there are always elements u(x) for which [w(a;),A] 7^ 0, 
so [G(u), G(X)} 7^ 0. Thus in nonzero field, only the components of Gauss' law in the 
Cartan sub-algebra are set to zero, and not the remaining components. 

This feature is generic to any effective potential which is a function of the Wilson 
line. Consider, for example, the potential for the real part of the trace of the Wilson 
line. In euclidean spacetime, 

Z tr (J) = [ VAJX,T) e S-jf<PxRe(trL) ^ ^ 

using a source J coupled to the trace of the Wilson line. Using the above techniques 
this equals the partition function 

Z tr (J) = f VVl(x) £ (M x )\ e-^-- 7 /^ 6 ^ \A?(x)) . (19) 

Again, while the sum over states is gauge invariant, the individual states are not: the 
source dependent state 

\Ai(x), J) = J VQ(x) e- J I d3xIie{tra) \A?(x)) (20) 

is not gauge invariant: \Aj , J) 7^ \Ai, J) when J 7^ 0. 

This behavior is not entirely unexpected. After all, the vacuum expectation value 
of the trace of the Wilson line, {tr L{x)) = dln(Z(J))/dJ(x)\j=o, is the wave function 
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for an infinitely heavy test quark [I]. In terms of Gauss' law, a single Wilson line 
corresponds to the insertion of a point source of unit strength at that point Jj]]. A 
source for the (trace of the) Wilson line, spread out over all of space, then modifies 
Gauss' law everywhere. From ([H]), when J(x) ^ no component of Gauss' law 
vanishes, not even those in the Cartan sub-algebra. Amusingly, the right hand side of 
Gauss' law, which is proportional to the source J, is only real if J is purely imaginary. 

Partition functions similar to that in ( |H| ) have been studied before jlQfl . The 
difference is that these works all assume that only gauge invariant states, as in (|17|), 
contribute: 

ZinvW = £ (Mx)\ mv e-P n+2itr ^\A(x)) mv } (21) 

Ai{x) 

At zero field everything reduces to the standard free energy: Z(0) = Z inv (0). Since 
Zi nv contains no term in the measure for the Vandermonde determinant, at one loop 
order Z(X) and Z inv (X) differ by a product of Vandermonde determinants, as the 
spatially longitudinal gluons give l/Yix l^(^)| 2 - This leads one to suspect that per- 
haps to arbitrary loop order, Z(X) and Z inv (X) differ only by terms which vanish in 
dimensional regularization. 

Nevertheless, (|2~TD has one serious drawback: in euclidean spactime it does not cor- 
respond to the functional integral for a gauge invariant potential. Using the definition 
of the invariant state \Ai(x))i nv , 

z inv (\) = fvn^x) fvn 2 (x) ]T (Mx)\ e-^ H \Af f (x)) , n f ( x ) = n 1 (x)e iX a^(x) 

Ai(x) 

(22) 

This is similar to the partition function in (|§D, with the replacement of the gauge 
transformation Q by Qf. In (|8|) the background field A is related to a gauge invariant 
quantity, the spatial average over the phases of the eigenvalues of Q(x). In (|22|). 
however, A is not related to anything gauge invariant, since the eigenvalues of Qf = 
Vtie lX Q} 2 are not e iX . Computing the euclidean path integral which corresponds to 
(f22|) |l|] leads to an expansion about constant Aq = X/g, which is gauge dependent 
beyond leading order. 

Consider the consequences if, by dint of prejudice, one insists that only gauge 
invariant states should contribute to the partition functions of gauge theories. Then 
construct the effective potential for the trace of the Wilson line, as in (|T9|) , everywhere 
replacing the states \Ai) with gauge invariant states, \Aj) inv . But as \Af) inv = \Ai) inv , 
this partition function gives a trivial potential for the trace of the Wilson line, 
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Zinv(J) = Tlx fDQ(x)e~ J f d3xRe ( trn ^. this potential is independent of temperature 
or the field content, and as an infinite product, vanishes in dimensional regulariza- 
tion. This is untenable, manifestly in contradiction with numerical simulations of 
euclidean lattice gauge theories. 

Partition functions as in (|T3) do produce novel thermodynamic behavior. In 



ordinary partition functions each state makes a positive contribution to Z, which 
implies that the pressure is positive. In nonzero field, however, as the sum over states 
includes a factor of cos(2 tr(XQ)), states with nonzero charge can make a negative 
contribution to Z(X), and consequently the pressure can be negative. For instance, 
even for a single scalar field coupled to an abelian gauge field, there are regions of A 
in which the pressure is negative. 

Dixit and Ogilvie observed that in theories with dynamical fermions, while 
Z(N) transforms of the usual vacuum have greater free energy than the trivial vac- 
uum, they are often metastable. The cosmo logical implications of such Z(N) "bubbles" 
were developed by Ignatius, Kajantie, and Rummukainen 0. In ref. |4j], however, 
it was pointed out that the pressure of the metastable Z(N) state can be negative, 
depending upon matter content of the theory. This includes the physically interest- 
ing case of three colors and six massless flavors, as applies at temperatures above the 
restoration of the SU(2) x U(l) symmetry. While the pressure of W, Z, and Higgs 
bosons acts to make the total pressure positive, it is rather disquieting to find that 
the pressure of a subsystem is negative. On these, and other grounds, the authors of 
argue that the metastable Z(N) states are not thermodynamically accessible. 

We suggest the contrary. Surely a system with negative pressure is not thermo- 
dynamically stable; indeed, it will not even be metastable for long. But we propose 
that if a path integral is gauge invariant and well defined in euclidean spacetime, then 
the corresponding partition function, as in (Tj) and (|T9"D, do represent physically re- 



alizable systems. It's just that there's a richer variety of partition functions possible 
in gauge theories. 

We thank M. Creutz, S. Gavin, C. Korthals-Altes, B. Svetitsky, and L. Trueman 
for helpful discussions. This research was supported in part by the U.S. Department 
of Energy under contract DE-AC02-76CH00016. 
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